Orbital and Spin contributions to the ^f-tensors in metal 

nanoparticles 



A. Cehovin 

Division of Solid State Theory, Department of Physics, 
Lund University, SE-223 62 Lund, Sweden 

CM. Canali 

Department of Technology, Kalmar University, 
391 82 Kalmar, Sweden, and Division of Solid State Theory, 
Department of Physics, Lund University, SE-223 62 Lund, Sweden 

A.H. MacDonald 

Department of Physics, University of Texas at Austin, Austin TX 18112 

(Dated: February 2, 2008) 

Abstract 

We present a theoretical study of the mesoscopic fluctuations of g-tensors in a metal nanoparticle. 
The calculations were performed using a semi-realistic tight-binding model, which contains both 
spin and orbital contributions to the (7-tensors. The results depend on the product of the spin- 
orbit scattering time rgo and the mean- level spacing but are otherwise weakly affected by the 
specific shape of a generic nanoparticle. We find that the spin contribution to the g'-tensors agrees 
with Random Matrix Theory (RMT) predictions. On the other hand, in the strong spin-orbit 
coupling limit St^q/Ti 0, the orbital contribution depends crucially on the space character of the 
quasi-particle wavefunctions: it levels off at a small value for states of d character but is strongly 
enhanced for states of sp character. Our numerical results demonstrate that when orbital coupling 
to the field is included, RMT predictions overestimate the typical g-factor of orbitals that have 
dominant d-character. This finding points to a possible source of the puzzling discrepancy between 
theory and experiment. 
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I. INTRODUCTION 



The difference in energy between the eigenvalues for the two spin states of an isolated 
electron in an external magnetic field B is Qq^-qB, where /ib is the Bohr magneton and 

is the so-called (^-factor for an isolated electron. If we ignore quantum electrodynamical 
corrections, it follows from the one-particle Dirac equation that Qq = 2. In an infinite crystal 
in the absence of external magnetic fields, time-reversal symmetry dictates that each electron 
Bloch state, labeled by k, is doubly degenerate (Kramer's degeneracy), provided that the 
crystal has inversion symmetry 0. This degeneracy is lifted by a magnetic field BB, and the 
energy difference between the two states can be expressed by means of a symmetric, 3x3, 
"^-tensor", 

Ae^ = fiBB{B^g^Bf\ (1) 
The tensor is related to the total magnetic moment of the state k (see below) 

/iB(/U )k = /UB(iv)k + /UB(fi'o5')k • (2) 

Here ^ibL and fiBdoS are the orbital and spin magnetic moments. The brackets (. . .)k 
denote the expectation value in the k-state. In the absence of spin-orbit interactions, the 
orbital angular momentum is quenched, and only the spin magnetic moment contributes 
to the level splitting. In this case the Kramer doublet is composed of two (opposite) pure 
spin states and the tensor ^k is isotropic, (^k)jj = 45jj. It follows that the (^-factor for a 
magnetic field in the z direction, defined as the square root of the tensor element (^k)zz, is 
equal to the free electron go. The effect of spin-orbit coupling is fourfold: (i) because the 
states forming a Kramer doublet are no longer pure spin states, their average spin is less 
than 1/2, which tends to decrease the typical (^-factor; (ii) the orbital angular momentum 
is no longer quenched and the corresponding magnetic moment can contribute to the level 
splitting, either by decreasing or increasing g; (iii) the tensor structure of Q is non-trivial 
(i.e. the response to a magnetic field is anisotropic and does not in general lead a moment 
aligned with the field); (iv) ^k can vary strongly with k. In bulk metals direct measurement 
of electron (yf-factors can be obtained only via conduction electron spin resonance (CESR) 
experiments. CESR involves transitions between states of the electron continuum. The two 
quantities measured in these experiments are the position of a resonance line, that is the 
average ^f-factor of the conduction electrons, over the Fermi surface, g^v, and the linewidth, 
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corresponding to the spin relaxation time Tso- In case of weak spin-orbit interaction the 
(7-factor shift Ag = g^v — go can be evaluated P| perturbatively in the spin-orbit coupling 
constant ^, Ag ~ C,/W, where W is the band width. Since W >> ^, the effect of spin- 
orbit interaction on bulk (^-factors is expected to be small, even for heavy elements like 
gold. Indeed experimentally g^v in Al is essentially equal to go, while for Au g^.^ = 2.1. 
When evaluated perturbatively, the spin relaxation time rgn can be related to Ag and to the 
"resistivity" relaxation time through the Elliott relationjsf, = This relation shows 
that scattering off impurities, surfaces, and phonons, which determine Tr, affects indirectly 
Tso, although the mechanism responsible for the spin relaxation is the spin-orbit interaction, 
which is essentially atomic in character. 

In this paper we focus on the (/-factors for electrons confined inside a metal grain of 
nanometer size, where the quantum energy spectrum is discrete and its individual quasipar- 
ticle energy levels e„ can be directly observed at low temperatures. By replacing the Bloch 
index k with the discrete quasiparticle orbital index n, equations ((T)) and (j21) describe how 
a Kramer doublet splits in a magnetic field. In contrast to the bulk case, the effects on the 
spin-orbit interaction on g - summarized in (i)-(iv) above - is expected to be enhanced, since 
the relevant energy scale with which the spin-orbit coupling strength ^ should be compared is 
not W but the much smaller single-particle mean- level spacing 6. In fact quantum finite-size 
effects on (/-factors in ensembles of metal nanoparticles have been investigated in the past 



experimentally via CESR. (For a review see Re^O and references therein.) To the best of our 
knowledge though, these experiments are riddled with many puzzling features, which make 



their comparison with proposed theories p, |a] very difficult. Furthermore in none of these 
experiments is there sufficient detail to extract information concerning statistical distribu- 
tion of (/-factors (see point (iii) above ). On the other hand, two groups 3, 0, S 3, using 
tunneling spectroscopy in single-electron transistors, have recently succeeded in measuring 
(/-factors of individual quasiparticle levels of a single metal nanoparticle. Measured (/-factors 
ranged from 0.1 to 2, depending on material, grain size and doping; (/-factors displayed large 
level-to-level fluctuations and strong dependence on the orientation of the applied magnetic 
field. Clearly (/-factors in small metal grains measured by single-electron spectroscopy have 
little in common with bulk (/-factors measured by CESR. 

The statistical properties and mesoscopic fluctuations of (/-factors measured experimen- 
tally in Refs. Q, 0, l^ll^ are well described by theoretical distributions based on Random 



3 



Matrix Theory (RMT), derived by two independent groups [ll|, ll3|. There is however a 
longstanding puzzle in the comparison between this theory and experiment. The (/-tensor 
distributions obtained from RMT are normalized to the average (((7^))^^^, which has to be 
evaluated by independent arguments. In the regime of strong spin-orbit interaction, the 



two theoretical mo dels [l^ predict (g^) to have contributions from both spin and orbital 
magnetic moments 

{g"^) = -^^soS + , (3) 
nil K 

where / is the elastic mean free path, R is the size of the particle and a is a constant of order 
1^^. In the limit of strong spin-orbit scattering t^qS — > 0, the spin contribution vanishes 
and onl y th e orbital contribution a{l/R) ~ ((L^))^ survives. The nanoparticles studied in 
Refs. y,^^ are not disordered and therefore / ~ R. Thus Eq. © predicts that {g^) should 
never be much less than 1. In noble-metal nanoparticles, however, the measured values of 
{g^) are typically between 0.05 and 0.1. 

RMT is a phenomenological approach which assumes that the statistical properties of 
interest depend only on the symmetry of the Hamiltonian. This assumption, although 
appealing and reasonable, is usually difficult to justifies rigorously. This fact, together with 
the discrepancy between theorv[l2. Q] and experiment mentioned above, motivates the 
theoretical study of (7-factors in metal nanoparticles presented here. Our investigation is 
based on a semi-realistic microscopic tight-binding model that we solve numerically. 

From the comparison between RMT and our microscopic theory we are able to shed 
some light on the discrepancy between the experimental measurements of {g"^) and the value 



predicted in Refs. ll^ The results of our model are in agreement with most of the RMT 
conclusions as to the functional form of the (7-tensor distribution and its dependence on the 
spin-orbit coupling strength. In particular we find that (/-tensors are strongly anisotropic 
when the nanoparticle shape is not perfectly symmetric, and that this anisotropy originates 
from mesoscopic quantum fiuctuations of the nanoparticle wavef unctions. However, we are 
able to go beyond RMT, in that our more detailed analysis allows us to clarify the role 
played by the orbital motion and demonstrate that its contribution to the average (/-factor 
is very sensitive to the character of the quasiparticle wavefunctions. Based on these results 
we conclude that the small value of {g"^) measured experimentally can be understood if the 
orbital angular momentum of the tunneling states, in contrast to the RMT assumptions, is 
still partially quenched. 



The paper is organized as follows. In Sec. HT] we introduce and discuss our model, with 
particular emphasis on the coupling between orbital motion and external field. In Sec. IIIII 
we illustrate our numerical results, we compare them with RMT predictions, and we dis- 
cuss their implications for the interpretation of the tunneling spectroscopy experiments. A 
summary and concluding remarks are presented in Sec. IIVI 



II. MODEL 



We model the nano-particle as a truncated fee crystal lattice with A''^ atoms. The shape 
of the system can be arbitrarily varied to simulate the variability of realistic nanoparticles. 
A spd tight-binding-model is used with 18 orbitals at each atomic site, including the spin 
degrees of freedom. 

The Hamiltonian, 

'H = 'Hha.nA + "^^SO + "^Zee, (4) 

has been introduced in a study of the quasiparticle properties in ferromagnetic metal 
nanoparticles jl^ . Here we give only a brief description of the terms in Eq.(j3]). 
The first term -ffband is an orbital part, 

"^band = 51 X! X! ^ ili,ii2,s'^tlJ.i,s^j,IJ.2,s i^) 

involving the Slater-Koster parameters [1^, ty^,f_i2,s- Atomic sites are labeled by i, j, and 
f'^ couples up to second nearest-neighbors. The indices /ii,/i2 label the nine distinct atomic 
orbitals (one 4s, three 4p and five 3d). The spin degrees of freedom, labeled by the index s, 
double the number of orbitals at each site. 

The second term describes a spin-orbit interaction, essentially atomic in character 

^80 = -^ 51 J2 {i^f^^s\L- S\i,fi',s')cl^^,Ci^^,^^,, (6) 

i fj,,fi',s,s' 

refiecting the fact that relativistic effects are important only when the electron is close to 
the nucleus. In spite of the local nature of these interaction, the effect of Tigo on the spin- 
orbit relaxation time rgo, compared to non- disordered infinite systems, is strongly enhanced 
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by the destruction of crystal symmetry due to the nanoparticle surface. This is, to a certain 
extent, similar to the mechanism described by the Elliott relation. 

A quantitative measure of the relative strength of the spin-orbit interaction is given by 
the dimensionless parameter 

The spin-orbit scattering is strong if A >> 1 and weak if A << 1. In the limit of 
weak spin-orbit interaction rgo can be calculated perturbatively by Fermi golden rule. Here, 
however, we use a more pragmatic approach: we define t^q in terms of the average spin-orbit 
quasiparticle energy shift [ij] 

^rso^{\^n-el\), (8) 

where e„ and e° is the n-th eigenvalue with and without spin-orbit interaction respectively 
and the average (...) is performed over the spectrum of the nanoparticle {3]. 

Since rgo decreases weakly with particle size L, while 6 varies as L~^, the effective strength 
of the spin-orbit interaction decreases with decreasing particle size. In the experiments 
of Ref. y, [l^ containing up to several thousand atoms, A can be as large as 10 for Au 
nanoparticles. In our theoretical studies we are able to deal numerically with nanoparticles 
containing only up to a few hundred atoms, which would yield A << 1 for a realistic value 
of ~ 100 meV. In order to achieve larger A values we therefore artificially increase the 
spin-orbit coupling strength ^. For a disordered dot of a generic shape our numerical results 
depend on the value of A but only weakly on the separate values of ^ and S. 

The Zeeman part Tizee in Eq. (jH) is conveniently divided into two terms, 

Hzee = -/iBE^-{ E {^,^^,s\g,S + L\^,^^',s')cl^,c,^^,^,,}+nZ'''' (9) 

i fi,iJ,',s,s' 

The first term of Eq. Q is the usual atomic contribution already considered in our earlier 
paper arising from the space dependence of the vector potential on an atomic length 
scale. In addition to this, there is a second term, a non-local orbital contribution arising from 
the magnetic flux encompassed by closed loops describing the paths of an electron hopping 
from site to site. Only the latter term is accounted for in the RMT approach, whereas the 
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former term has a larger importance in many respects. We discuss this contribution in detail 
in the next section. 



A. Non-Local Orbital Contribution (NLOC) 



The NLOC originates from delocalized electrons in an external magnetic field. The 



coupling to an external magnetic fie^ 



duction of the Peierls phases 
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d in tight-binding theory is accomplished by the intro- 



exp{iQij) = exp 



— / A{x) ■ dx 

hcJR, 



exp 



he 



(10) 



modifying the Slater-Koster parameters so that, 



Kl^"" = EE E (exp(^e,,) - i)t;;^,,„,c:,^,,c,,,,„ 



•^"^ c- fill 

i,j s /^i,/i2 

The non-trivial approximation of Eq. ()10|) . replacing the contour integral with a simple 
line integral has been justified by Ismail-Beigi et al. [l9|. 
Choosing the vector potential in the symmetric gauge, 



A = --fx5, (12) 

the Peierls phases can be conveniently rewritten in a form suitable for perturbation theory 
purposes. 



e 



27r 



[{R, - Ri) X {Ri + Rj)] ■ B. 



The phase factor of Eq. (fTUI) is expanded to first order in 



exp(ieij) = 1 - i^iBL^^^^ ■ B + 0{B 



where L^^oc jg defined as 



(13) 



(14) 



2ti 



-[{Rj-Ri) X {R, + R,)]. 



(15) 



4/iB0O 

This expansion allows one to write Eq. Q to first order in (small) B in the familiar form 



Hzee ~ - /iB fl- B = (16) 

-f^^Y.{ J2 {hf^,s\gsS + L\i,fi',s')cl^^,Ciy^^,} ■ B (17) 

i fi,fi',s,s' 
i,j fi,fi',s,s' 

We have explicitly checked that the Zeeman splitting due to the non-local interaction 
is gauge invariant in our numerical calculations, although the definition presented above is 
obviously for a particular gauge choice. 



B. Wannier states 

The basis set, |j, /i,s), is to be interpreted as a set of Wannier orbitals, extracted j3| 
from bulk properties. The Slater-Koster parameters, t^-J^^2,s' Hamiltonian matrix ele- 
ments of |j, /i, s). Application of t^fll^f^^^s ^ nano-particle rehes on their transferability, the 
assumption of local electronic environmental insensitivity to the boundaries of the system. 
The success of this procedure depends on the nature of the states involved . In noble metals 
the electronic density of the sp-bands is delocalized, and very different from a simple sum 
of atomic densities. The corresponding Wannier orbitals lose their atomic-like character in 
this limit, and the transferability is expected to be less robust, since the boundary condi- 
tions play an increasingly important role. The orbitals, \j,fi,s), of sp-symmetry will thus 
have large hopping elements, a property expected to remain the same in any transferable 
tight-binding treatment of these states. The (i-orbitals originate from bands of a different 
nature, with rather localized electronic densities around the atomic sites. The corresponding 
Wannier orbitals remain atomic-like and hopping parameters are small. 

From these arguments it is immediately clear that the contribution from 'hO^l'^^ becomes 
increasingly important in the limit of extended Wannier states with large hopping parame- 
ters. 

The eigenstates of this model are in general a mixture of all basis orbitals, | j, s) . An 
approximate separation of the eigenspectrum into bands of specific symmetries can still be 
made. For this purpose it is useful to define the projection operators, Ps, Pp and Pd- 
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Data 


Material 


<gi> 


< 52 > 


< 53 > 


A 


Exp. 


Cu 


0.9 ±0.3 


1.2 ±0.2 


1.6 ±0.3 


1.1 ±0.1 


RMT 




1.0 


1.1 


1.6 


1.1 


Here 


Au 


1.1 


1.3 


1.7 


1.1 


Exp. 


Cu 


0.4 ±0.2 


0.8 ±0.2 


1.3 ±0.3 


1.8 ±0.1 


RMT 




0.5 


0.8 


1.3 


1.8 


Here 


Au 


0.7 


0.9 


1.4 


1.8 



TABLE I: Calculated mean values of gi compared to random matrix theory predictions and ex- 
periments. The ^(-factors of our model are based on d-states only. The experimental and random 
matrix theory data are from Petta and Ralph [lo| . 



i S fig 

Pp = mm K/^P'S)(^>/^P'S| ' (20) 

Pd = ^^^\h^^d,s){i,|J,d,s\ , (21) 

with 

Ps + Pp + Pd = l. (22) 

These operators project onto the subspaces of s, p or (i-symmetry. As an example, in 
Fig. ^ we plot the expectation values of these projection operator in the eigenstates \n) 
of the Hamiltonian of Eq. (0}. The calculations are done for a 143 gold atom nanoparticle 
with hemispherical shape, from which 5 atoms have been removed to break the rotational 
symmetry. It is seen that states in the middle of the occupied band have well-defined 
(i-character, with the exception of some states around n = 700, which have mixed spd- 
components and low-energy states, which have mainly s-character. High-energy states far 
above the Fermi level have dominant p-character. States around the Fermi level are rather 
strongly mixed with the ci-component decreasing sharply with eigenvalue number. 
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n 



FIG. 1: Matrix elements of according to Eq. (|21() . The eigenstate |n > is for the Hamiltonian of 
a 143-atom hemispherical Au nanoparticle, with 5 atoms removed. The dashed vertical line marks 
the Fermi level. 

C. g-tensor 

In the absence of an external magnetic field, a nano-particle will have a doubly degenerate 
eigenspectrum, e„, reflecting the formation of Kramer pairs in a time- reversal symmetric 
system. In the presence of spin-orbit interaction, the Zeeman splittings, e„ — e„ ± (5e„, 
show anisotropic dependence on the direction of the applied field, B. For small B this is 
conveniently described in terms of the (7-tensor j^. 



By using standard degenerate perturbation theory in Hzee, the Gn matrices can be related 
to the matrix elements of the (dimensionless) magnetic moment operator fl, 




(23) 
(24) 



g, 



n 




(25) 



10 



{Gn)3,j = 2(V'„|/i,#„) . (26) 

Here \ipn) and {Tipn) are the time-reversed pair of eigenstates of 7i in the absence of B, 
corresponding to the eigenvalue e^. The diagonahzation of Qn, or in other words a suitable 
choice of coordinate system for the external field B, allows one to write 6en{B) in terms of 
the 3 eigenvalues g'^, j = 1, 2, 3, of Q 



S^niB) = ^^glBl + glBl + glBl (27) 

where Bj, j = 1, 2, 3, are the components of the magnetic field along the "principal axis" , 
61,62,63, defined by the three normalized eigenstates of Q. We refer to gj,j = 1,2,3 as 



the principal (^-factors ll|. The Zeeman splittings, 6en{B), can also be obtained from direct 
numerical diagonahzation of the Hamiltonian for a given direction of B. This approach 
is of course time demanding, but served as a check of the correctness of the perturbative 
approach, since the two methods should give identical results for small B. 

In the limit of zero spin-orbit interactions the three principal ^f-factors are identically 
equal to 2. In this case the Kramer pairs are pure spin states and the angular momentum is 
quenched. With the introduction of spin-orbit interaction spin characters mix, and non-zero 
orbital contributions {2^ to gi are expected. At weak spin orbit strength the spin character 
mixing is weak and the orbital contribution to g^ is small. For large A the spin character 
mixing is strong, while orbital contributions are expected to increase. 



III. RESULTS AND DISCUSSION 



In this section we present results of numerical calculations of (7-tensors for Au- 
nanoparticles containing 143 atoms. Similar results, not shown here, were obtained for larger 
and smaller sizes. We first consider the distribution of principal-axis directions. Fig.(j2K) 
shows calculated principal axis directions for a perfect hemisphere. The direction of is 
strongly dependent on the symmetry of the nanoparticle, with one axis always along the nor- 
mal of the base plane. Due to the crystal fee symmetry and the hemispherical truncation, the 
two remaining (7-factors will be degenerate and two principal axes can be arbitrarily chosen 
within the base plane. This is shown in Fig.(j2t), where we plot the direction of the principal 
axis directions on the unit sphere for all the quasiparticle orbitals of the nanoparticle. The 
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FIG. 2: Principal axis directions on the unit sphere for a 143 atom Au hemisphere in (a). In (b) 5 
atoms have been removed from the hemisphere. Each point represents one principal axis direction. 

situation changes dramatically with a small distortion of the shape obtained by removing five 
atoms from the hemispherical truncation. This slight breaking of the spherical symmetry 
produces a completely random distribution of principal axis, no longer preferentially aligned 
along crystalline or shape directions, as clearly seen in Fig. (j23). A random distribution of 
the principal axis directions was originally predicted by RMT ^ . Recent experiments Q| 
measuring principal gf-factors find random spatial orientation of the corresponding axis di- 
rections, Cj. The measurements were performed on Cu nanoparticles with approximately 
hemispherical shape. The results demonstrates that small irregularities in the boundaries 
in an otherwise ordered sample are enough to produce completely randomized directions, 
exactly as happens in our model. 

Once the perfect axial symmetry of the hemispherical nanoparticle is broken, for example 
by removing a few atoms, we find that the numerical results depend weakly on the particle 
shape: nanoparticles with shapes that are even less symmetric display qualitatively similar 
distributions of principal ^^-factors. It turns out that the dimensionless spin-orbit strength A, 
defined in Eqs. ©-(jHl), is the crucial parameter that controls the distributions. A meaningful 
evaluation of A(^) for a given value of ^ is however made complicated by the fact that states of 
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FIG. 3: A^(i^) in (a) shows the spin-orbit strength for a 143-atom Au nanoparticle with 5 atoms 
removed. XsiC) originates from a nanoparticle with a strongly disordered shape, but the same Na 
as A. The states are divided into d-states, < n|P^|n >> 0.7, and mixed states, < n\P(i\n >< 0.7 . 
(b) shows the differences in A between the two systems in percent. 

different orbital character, as determined by the projection operators discussed above, yield 
very different values of A for the same ^. Therefore we divide the eigenstates into two groups 
of states, which we refer to d- and mixed states. An eigenstate is operationally considered 
a (i-state if {n\Pd\n) > 0.7; it is a mixed state otherwise. We will use this distinction 
extensively below. For a given value of ^, the parameter A is then calculated separately for 
these two group of states. This procedure turns out to be very useful in the interpretation 
of our numerical results. A first example of this is shown in Fig. Q. In Fig. (jSt^) we plot 
A(^) as a function of ^ for these two groups of states for a given nanoparticle. In both cases 
A increases approximately linearly with ^, but the slope is much larger for the rf-states. In 
Fig. (jnb) we plot relative differences of the A for two nanoparticles with the same number 
of atoms {Na = 143) but very different shapes, again separating d- and mixed states. Here 
Xa refers to a nanoparticle with hemispherical shape with 5 atoms removed; Xb refers to a 
nanoparticle with a strongly disordered shape. It is seen that A for the d states is essentially 
the same for both nanoparticles for all values of ^. On the other hand mixed states give 
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FIG. 4: Mean square principal (^-factors of a 143 atom Au hemisphere with 5 atoms removed. The 
states are divided into d-states, < n|P^[n >> 0.7 , and mixed states, < n\P^n >< 0.7 . The 
full Hamiltonian is used in (a) while the NLOC to the Zeeman term is lacking in (b). The fitting 
function is given in Eq. H29() 

a value of A that is more sensitive on the shape of the nanoparticle, although the relative 
difference is always less than 10%. It follows that the two nanoparticles will have very similar 
(yf-factor distributions, with a stronger shape dependence of the ^f-factor of mixed character 
states. 

We now come to the discussion of the main results of our paper and to the comparison 
with RMT predictions. A very useful quantity for this comparison is the average of the sum 
of the squares of the principal (yf-factors 

1 



(28) 



In Fig. 0] we plot ((7^) as a function of A for a Au hemispherical nanoparticle with 5 
atoms removed, separating again the contribution from d- and mixed states. The results 
in Fig. m^a) are for the case in which the spin term and both orbital terms are included 
in the magnetic moment of Eq. (j2I). As expected, ((7^) ^ 4 in the limit A ^ of zero 
spin-orbit coupling and decreases monotonically with increasing A. In the limit of strong 
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spin-orbit coupling, A >> 1, {g^) tends to a constant value, which is relatively small (~ 0.3) 
for (i-states but is of order 2 for mixed states. This saturation value in the strong spin-orbit 
scattering limit comes almost entirely from the non-local orbital contribution (NLOC) to 
the magnetic moment. Indeed Fig. El^b) shows that that the saturation value is negligible 
when the NLOC is not included; in this case the {g"^) vs. A for d- and mixed states fall on 
the same curve which goes to zero like 1/A^ for large A. Note that the large saturation value 
for mixed states, when the NLOC is present, is consistent with the fact that rf-states are 
linear combination of atomic orbitals with small hopping parameters, whereas mixed states 
tend to be more free-electron like. 

We found that the calculated ((7^) (A) could always be fit to 



(/)(A) = V^ + 6^' (29) 



4-z/(A) 

TTa^ 

in the entire range of A. Here r/ is a constant independent of A, whereas z/(A) is a function 
weakly dependent on A such that v{\ = 0) = Qrf and z/(A) — ^ 1 for A >> 1. In practice 
a reasonable fit (shown in the figure) is obtained by fixing v = Qyf. Remembering the 
relationship between A and t^q given in Eq. ((7j), we can see that Eq. is perfectly 
consistent with the RMT prediction 12, in the strong spin-orbit scattering limit given 
in Eq. (jHl). The first term in Eq. comes essentially from the spin contribution to the 
magnetic moment. The second term originates from orbital contributions and it has been 
written in this specific way to make contact with the the RMT phenomenological parameter 
7] that describes its coupling of the magnetic field to the orbital angular momentum^^. The 
RMT parameter rf can be estimated fo r sp ecific physical systems by computing {n\Ll\n) 



in the absence of spin-orbit coupling |. For a ballistic sphere with diffusive boundary 



;d tor sp 

conditions one gets rf = l/8[l2,ll3- Thus for no disorder RMT predicts for {g'^) a saturation 
value ~ 1, which according to our model would correspond to the case of mixed states. 

We can further analyze the relationship between RMT and our microscopic model by 
comparing the distribution of the three individual principal (7-factors in different regimes of 
the spin-orbit coupling strength. In Fig. Eland Fig. [HI we plot the distribution of gi, g2, and 
(73 for d and mixed states respectively. Translating into RMT language. Fig. and Fig. [HI 
correspond to small and large orbital contribution respectively. The sub-cases (a), (b), (c) in 
each figure correspond to the regimes of weak, intermediate and strong spin-orbit scatte ring 



respectively. The values of A are chosen to compare with the cases considered in Refs. [ll , 

15 
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For weak spin-orbit interaction, we find that typically (71 ~ (72 < 2 and (73 ^ 2 when the 
spin contribution to the g tensor dominates (i.e. << A as in Figl^l^a) ). On the other 



IJ 



hand, when the orbital contribution dominates the ^f-tensor (r^ >> A as in Fig Efa) ' 
(72 ~ 2, (73 > 2. Both cases are in remarkable agreement with RMT predictions (ill. 

The trends of our numerical distributions for the cases of intermediate and strong spin- 
orbit scattering also agree with the RMT scenario Jll|. In particular in the strong spin-orbit 
scattering regime and for weak orbital contribution (see Fig. El^c) ), all three principal g- 
factors are peaked at values smaller than 0.5. In this case, typical values of g are always 
much smaller than the free electron value go = 2. 

We conclude this section by attempting a comparison between the results of our model 
and the experimental results of Ref. [lOl In Table H] we have summarized: (i) the experi- 
mental measurements of the average values of gi, i = 1,2, 3, for two Cu nanoparticles-raws 
labelled by "Exp"; (ii) the RMT results, including spin contribution only-raws labelled by 
"RMT" ; and (iii) the results of our calculations for rf-states, including both spin and orbital 
contributions -raws 

following. First the value of A corresponding to a particular nanoparticle is obtained^^ by 
matching the RMT value of {g'^) with the same quantity measured experimentally. Given 
this A, the average values of the principal (7-factors_(Qi), (^'2), and {gs) can then be calcu- 



lOl and the results are reported in Table 



lated theoretically by RMT. This was done in Ref. 
I. We can also compute the averages of gi within our microscopic model by choosing the 
parameter ^ for a given nanoparticle so that the value of A obtained from Eqs. [3IH1 is equal 
to the value extracted from the experiment using RMT. The averages predicted by the two 
theories (RMT and our macroscopic model) can then be compared with each other and with 
the values measured experimentally. 

From the table one can see that RMT - with orbital contributions neglected - predicts 
average values of g^ in good agreement with the experiments for two different A's correspond- 
ing to two different nanoparticles. The message that we want to convey here is that the 
average g values of rf-states predicted by our microscopic model, at the same nominal values 
of A defined through Eqs. [7||H1 are also in good a gree ment with RMT and the experiments, 
even when the orbital contributions are included 2^. Thus our resuls suggests a scenario, 
namely (/-factor distributions of states of mainly d character, which, if realized, would solve 
the puzzle posed by the comparison between RMT and experiment. Note that RMT is based 
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on a picture of single-particle wavefunctions essentially ergodic in space and it is therefore 
unable to capture the nature of wavefunctions more localized around atomic cores or defects, 
which instead are naturally described within our tight binding model. 

The problem with the scenario proposed here is that the Fermi level for noble-metal 
nanoparticles lies near states of both d-like and mixed-like character, (see Fig. Q ) and 
therefore the latter ones most likely play a very relevant role in the tunneling experiments. 
However if we use the procedure described above and compute (/-factor averages including 
mixed states, we find values that are too large in comparison with experiment. 

As possible solutions, we consider the following: (i) electrons are, in fact, tunneling into 
states of pure d character, (ii) Electrons are tunneling into states of mixed characters, but 
the nature of these delocalized states is strongly affected by the surface and therefore their 
characteristics, including their orbital angular momentum, are profoundly modified. This is 
an effect which is clearly not included in our model, and most likely it is more important for 
mixed states than for (i-states. (iii) electron-electron interaction beyond the simple mean- 
field approximation incorporated in the SK parameters could also modify the electronic 
states and partially quench their orbital angular momentum. This effect is also not included 
in our model. Note that recent studies 21] have demonstrated a sharp increase of electron- 
electron interaction due to surface induced reduction of screening. At the moment we find 
the occurrence of explanation (i) problematic: apart from the fact that in our calculations all 
states of pure d-character are below the Fermi level and therefore not available for tunneling, 
states of s character have much larger tunneling probabilities than (i-states, when (as in the 
present case) the barriers separating the grain from the electrodes are made of aluminum 
oxide. Explanations (ii) or (iii) or a combination of the two seem more compelling and 
alluring. However pursuing these lines requires theoretical modeling beyond the scopes of 
the present work. 

IV. CONCLUSIONS 

In summary, we have presented a theoretical study of the statistical properties of g- 
tensors for individual quasiparticle energy levels of metal nanoparticles, based on numerical 
calculations for a simplified but realistic model that treats spin-orbit interactions micro- 
scopically. Our theory of the (7-tensors includes both spin and orbital contributions to the 
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FIG. 5: g-factor distributions of d-states for three different values of A. The system is a 143 atom 
Au hemisphere with 5 atoms removed. 
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FIG. 6: ^(-factor distributions of mixed states for three different values of A. The system is a 143 
atom Au hemisphere with 5 atoms removed. 
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magnetic moment. We have shown that even small deviations from a perfectly symmetric 
nanoparticle shape cause random fluctuations in the quasiparticle wavef unctions, which, in 
turn, are the source of strong anisotropics of the g^-tensors and strong level-to-level fluctu- 
ations both in the principal (yf-factor values and in the directions of their principal axis. A 
dimensionless parameter measuring the strength of the spin-coupling controls the (yf-tensor 
distributions for nanoparticles of generic shape, in excellent agreement with the prediction 
of random matrix theory. Our work sheds light on the relative importance of the spin vs 
orbital contributions to the g^-factors and the strong dependence of the latter on the orbital 
character of the wave-functions. The presence near the Fermi energy of states of both d-like 
and sp-like character is responsible for aspects of the (yf-factor physics that arc not cap- 
tured by Random Matrix Theory. The small values of ^f-factors measured experimentally 
suggest that the orbital angular momentum of the tunneling states near the Fermi energy 
is most likely still partially quenched, even in the presence of strong spin-orbit scattering. 
Our calculations demonstrate that angular momentum quenching cannot be due simply to 
the strong d-orbital hybridization of states near the Fermi energy. Non-trivial changes in 
electronic structure might be responsible, perhaps due to irregularly shaped nanoparticle 
boundaries that produce orbitals more localized at the surface than those in the model that 
we have studied. Enhanced correlation effects near the surface, could also play a role. If 
surface imperfections are the main source of orbital localization, large g^-factors should be 
observable in nanoparticles with very regularly shaped boundaries in which all orbitals are 
extended across the entire particle. 
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